P‐wave transmission across fractures with nonlinear deformational behaviour     

X. B. Zhao  J. Zhao  J. G. Cai 《国际地质力学数值与分析法杂志》2006,30(11):1097-1112

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

Source signature and elastic waves in a half‐space under a momentary shear line impulse     

Moche Ziv 《国际地质力学数值与分析法杂志》2003,27(3):233-258

The transient deformation of an elastic half‐space under a line‐concentrated impulsive vector shear load applied momentarily is disclosed in this paper. While in an earlier work, the author gave an analytical–numerical method for the solution to this transient boundary‐value problem, here, the resultant response of the half‐space is presented and interpreted. In particular, a probe is set up for the kinematics of the source signature and wave fronts, both explicitly revealed in the strained half‐space by the solution method. The source signature is the imprint of the spatiotemporal configuration of the excitation source in the resultant response. Fourteen wave fronts exist behind the precursor shear wave S: four concentric cylindrical, eight plane, and two relativistic cylindrical initiated at propagating centres that are located on the stationary boundaries of the solution domain. A snapshot of the stressed half‐space reveals that none of the 14 wave fronts fully extend laterally. Instead, each is enclosed within point bounds. These wave arresting points and the two propagating centres of the relativistic waves constitute the source signature. The obtained 14 wave fronts are further combined into 11 disparate wave fronts that are grouped into four categories: an axis of symmetry wave—so named here by reason of being a wave front that is contiguous to the axis of symmetry, three body waves, five surface waves and two inhibitor waves—so named here by reason that beyond them the material motion dies out. Of the three body waves, the first is an unloading shear wave, the second is a diffracted wave and the third is a reflected longitudinal two‐branch wave. Of the two inhibitor waves, the first is a two‐joint relativistic wave, while the second is a two‐branch wave. The wave system, however, is not the same for all the dependent variables; a wave front that appears in the behaviour of one dependent variable may not exist in the behaviour of another. It is evident from this work that Saint–Venant's principle for wave propagation problems cannot be formulated. Therefore, the above results are valid for the particular proposed model for the momentary line‐concentrated shear load. The formulation of the source signature, the wave system, and their role in the half‐space transient deformation are presented here. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

The response of an elastic half‐space under a momentary shear line impulse     

Moche Ziv 《国际地质力学数值与分析法杂志》2003,27(3):207-232

The response of an ideal elastic half‐space to a line‐concentrated impulsive vector shear force applied momentarily is obtained by an analytical–numerical computational method based on the theory of characteristics in conjunction with kinematical relations derived across surfaces of strong discontinuities. The shear force is concentrated along an infinite line, drawn on the surface of the half‐space, while being normal to that line as well as to the axis of symmetry of the half‐space. An exact loading model is introduced and built into the computational method for this shear force. With this model, a compatibility exists among the prescribed applied force, the geometric decay of the shear stress component at the precursor shear wave, and the boundary conditions of the half‐space; in this sense, the source configuration is exact. For the transient boundary‐value problem described above, a wave characteristics formulation is presented, where its differential equations are extended to allow for strong discontinuities which occur in the material motion of the half‐space. A numerical integration of these extended differential equations is then carried out in a three‐dimensional spatiotemporal wavegrid formed by the Cartesian bicharacteristic curves of the wave characteristics formulation. This work is devoted to the construction of the computational method and to the concepts involved therein, whereas the interpretation of the resultant transient deformation of the half‐space is presented in a subsequent paper. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

地形作用下有限振幅Rossby波的稳定性分析     

周定文  彭晓林  罗德海 《成都信息工程学院学报》1990,(Z1)

本文用微分方程定性理论方法分析了地形对有限振幅Rossby波的稳定性影响。结果表明:地形北坡有利于波动出现不稳定,南坡使波动稳定;西坡使导式波动不稳定但使曳式波动稳定,东坡使导式波动稳定但使曳式波动不稳定。  相似文献   

Numerical prediction of blast‐induced stress wave from large‐scale underground explosion     

Chengqing Wu  Yong Lu  Hong Hao 《国际地质力学数值与分析法杂志》2004,28(1):93-109

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

Propagation of Extremely Low Frequency Waves Through the Ionosphere     

Singh  Ashok K.  Narayan  D.  Singh  R. P. 《Earth, Moon, and Planets》2002,91(3):161-179

The propagation features of extremely low frequency electromagnetic waves through the multicomponent ionospheric plasma are studied. It is shown that at relatively lower frequencies refractive index for right hand mode is higher than the left-hand mode, which is reversed at higher frequencies. The thermal temperature of plasma particle causes decrease in phase and group velocities of both right and left-hand modes. The crossover frequencies for different plasma models are computed and variation with ion concentration and thermal velocity is studied. Explicit expression for group velocity and travel time has been derived and studied numerically. Finally, we have presented simulation of the ion whistler spectrograms for Hydrogen, Helium and Oxygen ions present in the ionospheric plasma. The results are compared with the experimentally detected hydrogen and helium ion whistlers. The importance of the present study in the exploration of ionospheric plasma is illustrated.  相似文献   

沿海渔业天气警报服务系统     

王雪森  林朝昌 《成都信息工程学院学报》1993,(1)

本文简要地介绍了建立沿海渔业天气警报服务系统的必要性,系统的设计特点,现场试验情况以及今后建台布网的设想。特别是把短波单边带通信,地面波传播方式和双音报警技术引入天气警报服务系统,开辟了短波单边带气象通应用的新领域。  相似文献   

非线性Rossby波所满足的三阶Zakharov方程与低频振荡     

罗德海 《成都信息工程学院学报》1993,(4)

本文在基本气流具有水平切变的情况下,利用摄动法导出了非线性Rossby波所满足的三阶Zakharov方程,然后,考虑了基流具有弱切变的情况,通过使用三阶Zakharov方程研究了Rossby波列的第一类不稳定性问题。结果表明:通过非线性作用,大气中的Rossby波列可产生调制不稳定。同时,本文对这种不稳定的区域,增长率和周期进行了详细的计算,并讨论了波振幅、波数、纬度和基流切变对它们的影响,指出Rossby波列的调制不稳定可以激发30～60天的低频振荡。  相似文献   

Precursory arrivals to PP     

Mamoru Kato  Kazuro Hirahara 《Geophysical Journal International》1991,106(3):551-557

  相似文献   

Accuracy of kinematic wave and diffusion wave approximations for space-independent flows on infiltrating surfaces with lateral inflow neglected in the momentum equation     

V. P. Singh 《水文研究》1995,9(7):783-796

Error equations for the kinematic wave and diffusion wave approximations with lateral inflow neglected in the momentum equation are derived under simplified conditions for space-independent flows. These equations specify error as a function of time in the flow hydrograph. The kinematic wave, diffusion wave and dynamic wave solutions are parameterized through a dimensionless parameter γ which is dependent on the initial conditions. This parameter reflects the effect of initial flow depth, channel-bed slope, lateral inflow, infiltration and channel roughness when the initial condition is non-vanishing; it reflects the effect of bed slope, channel roughness and acceleration due to gravity when the initial condition is vanishing. The error equations are found to be the Riccati equation. The structure of the error equations in the case when the momentum equation neglects lateral inflow is different from that when the lateral inflow is included.  相似文献